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The ground state of the one-dimensional t-J model coupled with phonons in the adiabatic limit 
is numerically investigated by use of the Lanczos technique at quarter filling. Due to the inter- 
play between the electron-electron Coulomb repulsion and electron-phonon interaction, this model 
shows a sequence of lattice distortions leading to the formation of charge-density-waves and bond- 
order-waves. Moderate electron-electron and electron-lattice coupling may lead to coexistence of 
dimerization and tetramerization in the distortion pattern. Dimerization leads to the formation of 
an "antiferromagnetic" Mott insulator, while tetramerization gives rise to a spin-Peierls phase. By 
increasing the super-exchange coupling, antiferromagnetism is inhibited due to the change of the 
distortion periodicity. 

PACS numbers: 75.10.Jm 05.70.Jk 75.40.Cx 



I. INTRODUCTION 

In the last decades the physical and chemical proper- 
ties of quasi one-dimensional (ID) systems have attracted 
growing attention due to the occurrence, in the same 
sample, of several different phases, ranging from insulat- 
ing to superconducting, when a control parameter, like 
pressure, is variedi^. Inorganic compounds, as CuGeO^ 
or iVaV^Og^, show clear evidence in favor of the existence 
of spin-Peierls states with broken lattice symmetry. Even 
more complex phase diagrams are displayed by quasi ID 
organic materials: the DI—DCNQI compounds^, the 
Bechgaard salts series (TMTSF^A and its sulfur ana- 
log (Fabre salts series) (TMTTF) 2 X, where X denotes 
different monovalent anions (X = PF§, AsFq,CIOa, ...). 
These systems, which consist of stacks of weakly cou- 
pled chains, exhibit, at low temperature, several bro- 
ken symmetry phases like charge-density- wave (CDW), 
spin-density-wave (SDW) and spin-Peierls (SP) before 
becoming superconducting^. The physical mechanism 
underlying these strongly different orderings can be iden- 
tified as the competition between lattice distortions and 
electron dynamics, with a major role played by elec- 
tron interactions. Unfortunately, interaction effects in 
low dimensional electron systems coupled to further de- 
grees of freedom, such as impurities or phonons, are not 
fully understood yet. Moreover, in order to describe the 
phase diagram of these materials the inclusion of inter- 
chain coupling besides interaction and phonons cannot be 
avoided. At the very least, inter-chain coupling is crucial 
for the stabilization of a genuine three dimensional long 
range order, like antiferromagnetism or superconductiv- 
ity However, we believe that, as a first step, a thorough 
understanding of the behavior of purely ID strongly cor- 
related electron systems coupled to adiabatic phonons is 
necessary. 

By now, the physics of one dimensional electron sys- 
tems in a rigid lattice is well understood^: spin isotropic 
models display a metallic phase, the celebrated Luttinger 



Liquid (LL), with gapless charge and spin excitations. 
Umklapp processes occurring at commensurate fillings 
drive the opening of a charge gap, while attractive in- 
teractions may lead to pairing and then to a spin gap. 
The presence of a gapless branch of excitations is inti- 
mately related to the long range power law behavior of 
correlations and to the divergence of selected susceptibili- 
ties in these models. Remarkably, in ID all these features 
are governed by a single parameter, denoted by Kfi. In 
the LL phase both CDW and SDW susceptibilities at 
wave vector q = 2fcp diverge for K p < 1 while the CDW 
susceptibility at q = Akp is singular for K p < 1/2. Even- 
tually, superconducting fluctuations become relevant for 
K p > 1. The K p parameter has been numerically eval- 
uated as a function of the coupling constants in several 
ID models^. If the Luttinger Liquid is weakly coupled to 
classical phonons, second order perturbation theory pre- 
dicts the occurrence of lattice distortions of modulation 
q whenever the associated CDW susceptibility diverges. 

In order to substantiate these perturbative results and 
to extend them to finite phonon coupling, the numeri- 
cal study of simple lattice models on finite chains may 
indeed be profitable. The repulsive extended Hubbard 
model has been extensively investigated^iiSiiiiiSii^ and 
the interplay between the 2kp and Akp instabilities has 
been discussed. In such a class of models, the antiferro- 
magnetic correlations are generated by virtual hopping 
processes and do not lead to an effective attraction be- 
tween electrons. As a consequence, the key parameter 
K p is restricted to the "repulsive" range K p < 1. By 
contrast, in the t-J model the magnetic coupling J is 
unrelated to the on site repulsion (U = oo) and K p ac- 
quires values in the extended domain 1/2 < K p < ooi4. 
An investigation of the ID t-J model coupled to adiabatic 
phonons would then allow to single out the role of antifer- 
romagnetic fluctuations in driving lattice distortions and 
stabilizing CDW, SDW or SP ground states. This is the 
program we intend to pursue in the following sections. 
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II. THE MODEL 

One of the simplest models able to describe the low en- 
ergy properties of strongly correlated electron systems is 
the t-J model. Here, contrary to the Hubbard model, the 
magnetic interaction is explicitly included in the hamil- 
tonian while the strong on site Coulomb repulsion acts 
via the single occupancy constraint. The phase dia- 
gram of the t-J model in a rigid ID lattice has been 
determined^ and analytic solutions are known in the 
J — ► limit (spinless-fermion case)i£ and at J/t = 2 
(super-symmetric case)i£*i£, besides at half filling, i.e. for 
electron density n = 1, where it reduces to the Heisenberg 
chain. In the latter case the model is known to be unsta- 
ble toward dimcrization (spin-Pcicrls instability) when 
coupled to an elastic latticed. As a result, a periodic 
modulation of the lattice at wave- vector q = 2k p = t sets 
in and the system displays a bond-order wave (BOW): 
adjacent spins pair in a singlet state and a gap in the 
spin spectrum opens up. When phonons are treated in 
the adiabatic limit, this transition has been predicted for 
arbitrarily weak spin-lattice coupling— 

At finite doping, i.e. for electron density n < 1, the 
pure ID t-J model enters the Luttinger Liquid phase 
which extends up to J/t < 2.3 (for n = 1/2) with growing 
superconducting fluctuations until, at larger J/t, phase 
separation sets in. Contrary to the purely repulsive Hub- 
bard model, the t-J model explicitly exhibits the attrac- 
tion mechanism between electrons mediated by antifer- 
romagnctic fluctuations which is likely to be present in 
real magnetic materials. 

In the following, we study ID t-J model coupled with 
classical phonons and we investigate by numerical diag- 
onalization technique the interplay between electron in- 
teractions and lattice distortions. We focus our attention 
on the competition, or sometimes the cooperative behav- 
ior, between charge ordering and lattice instabilities and 
we examine its consequences on the magnetic properties 
of the materials. The calculations will be performed at 
quarter filling (n = 1/2), which is the appropriate choice 
to mimic the behavior of the organic compounds^, and 
our results are summarized in Fig. 1 which shows a ten- 
tative zero temperature phase diagram of the model. 

The one-dimensional t-J model coupled with adiabat- 
ically phonons is defined by the hamiltonian 

H = - (1 - Si) (cl t(r Ci+x,a + h.c?) + 



(1) 



Where Si are spin- 1/2 operators at the site i, ct = 
cj a (1 — ni.-a) are electron Gutzwiller-projected cre- 
ation operators and m = J2a c l a c i,<y- We have also set 



the bare hopping integral t to unity thereby fixing the 
energy scale. The hamiltonian depends on the classical 
variables Si which identify the bond distortions defined 
as Si = (u;+i — Ui) where m is the displacement of ion i 
with respect to its equilibrium position. The last contri- 
bution represents the elastic deformation energy and Kb 
is the spring constant. The form of the electron-phonon 
coupling adopted here follows from a linearization of the 
expected dependence of the coupling constants on lattice 
distortions: for small displacements 5 of the ion sites, 
the hopping and exchange integrals vary with the dis- 
tance as a certain power— a which has been estimated 
in the range 6 < a < 14 



Ji,i+i (8) = J 



■ U i+ l 



J 1 
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where a is the bare lattice spacing. For simplicity we have 
taken the same a for both the hopping and the super- 
exchange integrals. We note that the Hamiltonian (1) 
is invariant under the rescaling a — > Xa* , Kb — ► X 2 K* Bl 
5 — > S*/X, where A is the rescaling factor. We can then fix 
one parameter between a and Kg2^ without affecting the 
physics of the model. In the hamiltonian JIJ we have set 
a = 1 . Clearly, the physical range where the hamiltonian 
||TJ may be used is restricted to \Si\ -C 1. 

This model is studied at quarter filling (n = 1/2) and 
zero magnetization by use of the Lanczos method on a 
finite ring 16 sites . Contrary to previous studies on the 
Hubbard model, we choose open shell boundary condi- 
tions, i.e. periodic boundary conditions when the number 
of lattice sites L is an integer multiple of 8, anti-periodic 
boundary conditions otherwise. This choice is motivated 
by the attempt to mimic the behavior of an infinite chain 
in a small ring. As already noticed, a key feature of the 
t-J model in the thermodynamic limit is the presence of 
gapless excitations which lead to a diverging CDW sus- 
ceptibility, and then to lattice distortions. In closed shell 
rings, this feature is lost because of the 0(1/ L) finite size 
gap in the kinetic energy spectrum. As a consequence, 
unrealistically low values of the elastic constant Kb are 
needed in order to stabilize lattice distortions— 

The ground state of the hamiltonian Q is numerically 
obtained by Lanczos technique at fixed {Si}. The lat- 
tice distortions are then updated in order to minimize 
the total energy E({Si}). This can be accomplished by 
solving the self-consistency equations which follow from 
the Hcllmann-Feynman theorem: 



dE 
dS~ 



y^(cf,<T C i-l-l,g + h - c ) + Jl \ S i ' &+! - jTHTli+l) 

' S i+ 1 ^ \ n 3 n J+^l = ( 2 ) 



The constraint J^i^i = which reflects the chosen ring 
geometry has been implemented by use of a Lagrange 
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multiplier. Here (...) is the ground-state expectation 
value. In order to find the bond length configuration 
which minimizes the total energy at given coupling con- 
stants (J, Kb) we iterate the set of equations © until 
convergence is reached. 



III. NUMERICAL RESULTS 

Here we present the numerical results on the t- J model 
coupled to adiabatic phonons at quarter filling and van- 
ishing magnetization for a 16 site ring. In this case 
fci? = mr/2 = 7r/4 and the two expected, competing pe- 
riodicities correspond to tetramcrization (q — 2k F ) and 
dimcrization (q = Akp). In the former case two electrons 
belong to the unit cell and then we expect a paramagnetic 
insulator behavior with both charge and spin gap. This 
regime precisely corresponds to the spin-Peierls phase ex- 
perimentally found in the Fabre salt series. Instead, in 
the presence of lattice dimerization, the model should be 
a paramagnetic metal, having a single electron per unit 
cell. However, the residual Coulomb repulsion drives the 
system toward a Mott insulator phase. Such a circum- 
stance can be explicitly checked in the J — * limit, where 
some exact statement can be made. 

In one dimension, when the kinetic term alone is 
present, the charge degrees of freedom arc described by a 
free spinlcss fermion gasi^ whose properties, even in the 
presence of (adiabatic) phonon coupling, can be explic- 
itly derived in the thermodynamic limit. Lattice dimer- 
ization is stabilized for arbitrarily large KgQ but such a 
symmetry breaking cannot be accompanied by a CDW 
of the same periodicity because reflection across a bond 
leaves the distortion pattern unaltered while interchang- 
ing the two sublattices. As a result, in the t-J model at 
n = 1/2 and J — > 0, the charge distribution remains uni- 
form along the chain while a charge gap opens due to the 
doubling of the unit cell. The spin degrees of freedom 
arrange as in the ground state of the Heisenberg chain 
and the spin correlations can be expressed in a factor- 



(S r ■ So) 



r+l 

Y,Psf U)Sh(j-I) 
i=2 



(3) 



where Pg F (J) is the probability of finding j particles in 
(0, r) with one particle in and another in r, evaluated 
in the ground state of the dimerized spinlcss Fermi gas. 
If we set N r = X)[=o n i' the probability Pg F (j) can be 
formally expressed as Pg F (j) = {n n r d (N r — j)} and 
ShU) = (Sj • &o)h is the spin-spin correlation function 
of the isotropic Heisenberg chain. Following Ref^£ the 
asymptotic decay of Eq. © in the presence of lattice 
dimcrization can be evaluated as 



(S r • S ) oc 



cos(2fcpr 



(4) 



much slower than in the standard t-J model at the same 
electron density. This asymptotic result valid for any 




FIG. 1: Phase diagram of the t-J model coupled to adiabatic 
phonons obtained by numerical diagonalization of a 16 
site chain. The boundaries of the uniform phase (U), the 
dimerized BOW (4ft j?), the tetramerized BCDW (2k F ) and 
the phase separated region are shown by full lines. The dashed 
line represents a crossover between the bond pattern shown 
in Figs. 2c (left) and 2d (right). 



lattice dimerization, can be understood from Eq. |J3J| 
the simpler limit of decoupled dimers, i.e for Si ~ (— 1) ! 
where 



in 

+i 







PsfU) 



fj,(r+l)/2 
5 3 ,(r+2)/2 



for r = 1 
for r odd ^ 1 
for r even 



(5) 



Remarkably, Eq. Q reproduces the same 1/r decay of 
spin correlations as in the usual Heisenberg chain. This 
form of the magnetic correlations, together with the pres- 
ence of a charge gap shows that the model is indeed in a 
Mott insulator phase. Such a result suggests that lattice 
dimcrization strongly favors antifcrromagnctism, which 
will be eventually stabilized by inter-chain coupling giv- 
ing rise to three dimensional antiferromagnetic ordering. 
Therefore, in the J — > limit, the t-J model coupled to 
adiabatic phonons reproduces one of the phases experi- 
mentally found in the quasi one dimensional Fabre salts. 
Now we want to investigate how the enhancement of an- 
tiferromagnetic coupling induced by the J term modifies 
this picture. 

The numerical results are summarized in Fig. 1 
where the phase diagram emerging from our analysis is 
sketched. No lattice or spin symmetry has been imposed 
and the Lanczos diagonalization has been performed in 
the full Hilbert space of the model. This strongly limits 
the feasible sizes to the range L < 16. In all cases the 
ground state we found shows definite lattice periodicity, 
the result does not depend on the starting bond config- 
uration or on the convergence requirements. Noticeable 



exceptions to this statement are the states in the phase 
separated region (PS in Fig. 1) at very large J, where the 
energy landscape presents several local minima. In order 
to check the stability of the numerical results, we per- 
formed a set of Lanczos diagonalizations in the reduced 
Hilbert space obtained by selecting a given unit cell (e.g. 
two sites or four sites). We always found consistency be- 
tween the numerical output in the full and in the reduced 
Hilbert space except in the PS region, where the ground 
state energy obtained in every subspace of definite trans- 
lational symmetry is larger than the ground state energy 
found in the full Hilbert space. This feature identifies 
the boundary of the PS region. We will not discuss this 
unphysical regime any further. 

The distortion patterns we found correspond to peri- 
odicities of growing wavelength moving from small J to 
larger J: 

- uniform (U in Fig. 1) 

- dimerized (4ftp in Fig. 1) 

- tetramerized (2k p in Fig. 1) 

A general feature of the phase diagram we obtained 
is the rather weak dependence of the phase boundaries 
on the spring constant Kg: this shows that the antifcr- 
romagnetic coupling J is indeed the driving force which 
selects the periodicity of lattice distortions. In the uni- 
form phase, Si = and the system is a metal according 
to the theory of Luttingcr Liquids. The numerical data 
predict a uniform phase also in the previously discussed 
J — > and large Kb limit where the charge degrees of 
freedom behave as spinless fermions which are known to 
undergo spontaneous dimerization when coupled to adi- 
abatic phonons. This shows that, at least in this region, 
strong finite size effects are present in our small cluster. 
In order to investigate this problem, we performed few 
Lanczos diagonalizations for the L = 20 site chain in 
the Hilbert sub-space defined by the eigenstates of the 
four site translation operator. This calculation allows 
to study the interplay between uniform solution, dimer- 
ization and tetramerization in this cluster. The phase 
boundary between the uniform and the dimerized phase 
shrinks from J ~ 1.2 to J ~ 0.8 at If # = 3 while the on- 
set of tetramerization hardly moves. Therefore we arc led 
to conclude that, as already discussed in previous numer- 
ical studies of the Hubbard modelii, the very existence 
of the uniform phase in the adiabatic limit is likely to 
be an artifact of the smallness of the cluster considered, 
while the other phase are much less affected by finite size 
problems. 

At larger J the dimerized Akp state prevails and ac- 
cording to the previous analysis the chain behaves as a 
Mott insulator. While the electron density remains uni- 
form, the lattice develops a bond ordered wave (BOW): a 
charge gap opens due to the Umklapp processes triggered 
by the doubling of the unit cell and the cell-cell spin cor- 
relations show an extremely slow decay characterized by 
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FIG. 2: Pictures of the competing density waves in the 
ground state, a) Uniform state (U), b) Ak F BOW and uniform 
charge, c) 2k F BCDW tetramerized phase (S-W-S-W), 
d) 2k F BCDW tetramerized phase (U-S-U-W). 

a diverging antiferromagnetic susceptibility. When three 
dimensional coupling is allowed for, a genuine antiferro- 
magnetic long range order sets in. 

Finally, by further increasing J, the lattice periodic- 
ity changes and the unit cell doubles again giving rise 
to lattice tetramerization. In this case, a weak 2kp = 
it/ 2 charge density wave develops on top of the lat- 
tice distortion giving rise to a bond/charge density wave 
(BCDW)±2iii: 

(rii) = g + A ^f cos + ^) 

Si = Bik P cos (m) + B 2 k F cos (J^ij (6) 

The phase shifts shown in Eq.© are those obtained by 
the numerical diagonalization and correspond to a CDW 
of the type — h + — , where + (— ) identifies a local den- 
sity higher (lower) than the average (n = 1/2). No sig- 
nificant Akp component in the CDW is detected in the 
BCDW, in agreement with previous investigations of the 
Hubbard model. Interestingly, the bond pattern depends 
on the relative magnitude of the B±k F and B2k F ampli- 
tudes: For Bik F > \B2k F \ the ordering is of the type 
S-W-S-W (Fig. 2c) while for \B Akp \ < B 2kp the 
sequence is U — S - U -W (Fig. 2d). Here, W, U and 
S respectively represent a weak (W) almost undistorted 
(U) and strong (S) bond. Although no real transition, 
but rather a crossover, is present between these two cases, 
we can identify a small region (shown in Fig. 1) where 
the bond ordering is of the first type, while in the re- 
maining region of stability of the tetramerized phase the 
bond lengths arrange as in Fig. 2d. The crossover be- 
tween these two different realizations of the 2kp phase 
may be understood by noting that the S — W — S — W 
BOW (Fig. 2c) appears as an intermediate regime be- 
tween the 4kp region and the stable U — S — U — W (Fig. 
2d) ordering, which indeed dominates in a large portion 
of the phase diagram. 

Notice that previous studies of the extended Hubbard 
modeliii^ found only a pattern of the first type, which is 
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FIG. 3: Energy landscape at J = 0.7 and K B = 1. The 
ground state energy, referred to the undistorted state, is plot- 
ted as a function of the difference A between the two weak 
bonds in the S — W — S — W' phase. The strong bond is fixed 
at the equilibrium value S = 0.15925 



often referred to as the dimerization of an already dimer- 
ized state. In such a regime, our calculations suggest 
the occurrence of a transition between a low tempera- 
ture tetramerized phase and dimerization at higher tem- 
perature. In order to investigate this possibility, which 
is indeed experimentally found in DI — DCNQI com- 
pounds, we evaluated the ground state energy of the 16 
site ring at given lattice distortion: We kept the am- 
plitude of the two strong bonds fixed while varying the 
difference A between the two weak bonds of the unit cell 
S-W-S-W. The result, shown in Fig. 3 for J = 0.7 
and Kb = 1, displays a double well structure character- 
ized by a small energy barrier between the tetramerized 
ground state at A ~ ±0.06 and the dimcrized configura- 
tion corresponding to A = 0. 

Whenever lattice tctramcrization prevails, two elec- 
trons belong to the unit cell and then the chain becomes a 
paramagnetic band insulator as in the spin-Peierls phase. 
We believe that, because of the presence of both charge 
and spin gap, this phase is stable toward 3D coupling. 

It is interesting to investigate the energy difference be- 
tween the electronic states corresponding to the low en- 
ergy lattice distortions which are solutions of Eq. 
Such a calculation may provide some useful information 
about the temperature range where an ordered phase is 
indeed stable. In Fig. 4 we plot the energy of different 
distortions as a function of J at fixed Kb = 3. More pre- 
cisely, we plot the energy difference between the lowest 
metastable state of a given periodicity (4k p or 2k p) and 
the uniform solution. 

The numerical data also suggest first order boundaries 
between the uniform and the dimcrized phase, as can be 
inferred by the sudden jump of the appropriate Fourier 
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FIG. 4: Energy gains between different distortions at fixed 
Kb = 3 in a L — 16 site ring. Open (full) squares refer to a 
four (two) site unit cell. 




FIG. 5: Amplitude of the 4fc_F (full dots) and 2kp (open dots) 
Fourier components of the bond distortions in the ground 
state of the model. Upper panel: Kb — 3, lower panel Kb — 1 



component of the lattice distortion {5i} shown in Fig. 5 
for a couple of values of Kb- Our results do not allow 
to draw any definite conclusion about the order of the 
transition between the dimerized and the tetramerized 
phase. 

Although the periodicities of the lattice distortions we 
found numerically precisely correspond to those expected 
on the basis of the instabilities of the homogeneous t-J 
model, some discrepancy on the sequence of distortion 
patterns does occur. As previously discussed, the t-J 
model is characterized by a K p parameter always larger 
than 1/2 in the whole range J > 0. This would imply 
that, at least in the perturbative regime of small distor- 
tions (i.e. for Kb — > oo), pure dimerization should never 
be observed. Instead, tetramcrization should be seen for 
< J < 2.3 which, according to Refii^, corresponds 
to K p < 1. While the upper limit for 2kp distortions 
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roughly agrees with the Lanczos diagonalization data, the 
numerical results clearly show that, for large Kg, dimer- 
ization and not tetramerization provides the most stable 
distortion pattern. This behavior can be understood by 
noting that the amplitude of the 2Uf singularity in the 
charge response of the uniform t-J model is strongly sup- 
pressed for J < 2, as already noted in Refsi 2 ^ 2 ^. This 
suggests that even if in the thermodynamic limit a 2kp 
singularity may eventually prevail in the whole phase di- 
agram, at small J a clear dimerization pattern with only 
weak tetramerization corrections will be visible. 



IV. CONCLUSIONS 

In summary, we determined the phase diagram of the 
ID t-J model adiabatically coupled to the lattice by use 
of Lanczos diagonalization in small clusters. We found, 
for increasing antifcrromagnctic coupling, a metallic uni- 
form phase, a BOW phase which displays the typical 
features of a Mott insulator and preludes to 3D anti- 
ferromagnetic (AF) order, a BCDW phase which realizes 
the spin-Peierls scenario. Therefore, by increasing the 
AF coupling J in the model the systems goes from anti- 
ferromagnetic to paramagnetic because of the change in 
distortion pattern determined by the charge degrees of 
freedom. This shows that in the t-J model charge and 
spin ordering compete due to the coupling to the lattice 



and strong super-exchange interaction may lead to a spin- 
Pcicrls ground state. Remarkably, the phases we found 
are all present in the typical zero temperature phase dia- 
gram of the quasi ID organic materials^. Moreover, these 
phases do indeed occur in the same sequence as we found 
when uniform pressure is applied to the material. The 
effect of pressure is mainly to increase the hopping ampli- 
tude (which we set as unit of energy) and then our results 
are compatible with the experiments provided the effec- 
tive intra-chain supcrexchange coupling J is only weakly 
dependent on the applied pressure. If this is the case, 
the t-J model may indeed represent a useful theoretical 
framework for the interpretation of the behavior of the 
quasi ID organic materials. 

The effects of charge ordering in this class of materi- 
als has been mainly investigated by using an extended 
Hubbard model coupled with the lattic e) 11 ' 12 . A phase 
diagram, of the ground state, obtained as a function of 
the electron-phonon couplings show the same periodici- 
ties that we found with the exception of the one in Fig. 
2d. Notice that the strong coupling limit of the extended 
Hubbard model is related to an effective t-J modcl 2 ^ by 
J = u _ v , which restricts the mapping to the J <C 1 
regime. 
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